Synchronization in coupled map lattices as an interface depinning 
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We study an SOS model whose dynamics is inspired by recent studies of the synchronization 
transition in coupled map lattices (CML). The synchronization of CML is thus related with a 
depinning of interface from a binding wall. Critical behaviour of our SOS model depends on a specific 
form of binding (i.e., transition rates of the dynamics). For an exponentially decaying binding the 
depinning belongs to the directed percolation universality class. Other types of depinning, including 
the one with a line of critical points, are observed for a power-law binding. 
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Recently synchronization of chaotic systems attracted 
a lot of interest T|. To a large extent, this interest is 
motivated by numerous experimental realizations of this 
phenomena including lasers, electronic circuits, or chem- 
ical reactions 2]. Synchronization acquires additional 
features in spatially extended systems, where it can be 
regarded as a certain nonequilibrium phase transition. 
There are increasing efforts to understand the proper- 
ties of this transition. Relatively well understood is the 
synchronization transition (ST) in certain cellular au- 
tomata. Since a synchronized state can be regarded as 
an absorbing state, for cellular automata, which are dis- 
crete systems, the phase transition, as expected, belongs 
to the Directed Percolation (DP) universality class Q- 
However, continuous systems, as e.g., coupled map lat- 
tices (CML) 0, need infinite time to reach a synchro- 
nized state and the relation with DP does not seem to 
hold. Indeed, Pikovsky and Kurths argued that for 
continuous systems this transition should belong to the 
so-called bounded Kardar-Parisi-Zheng (BKPZ) univer- 
sality class 0. Recently, precise numerical calculations 
confirmed their predictions but only for CML with some 
continuous maps 0- Surprisingly, the ST for discontin- 
uous 01 13 ) or continuous but sufficiently steep maps 
was found to belong to the DP universality class. It 
would be desirable to understand the critical behaviour 
of the ST in CML and the present paper might be a step 
in this direction. 

Let us briefly describe the setup which is used to 
study synchronization in CML ^7]. In the simplest, one- 
dimensional case, one takes a single-chain CML of size 
L that is composed of L diffusively coupled local maps 
f{ui) 4]. The maps are chaotic and act on continu- 
ous site variables Ui {i ~ 1 , . . . , L) , that are typically 
bounded {0 < Ui < 1). Then, one couples such a spatio- 
temporally chaotic system with its identical copy, which 
initially has a different set of site variables. It turns out 
that the evolution of such a coupled system depends on 
the coupling strength. For weak coupling the two CML's 
are desynchronized and essentially independent. How- 
ever, for sufficiently strong coupling the system gets syn- 
chronized and approaches a state where corresponding 
pairs of site variables in both copies take the same val- 
ues. To quantify synchronization one can introduce the 
synchronization error w{i,t) = \ui{i,t) — U2{i,t)\ where 



a lower index denotes a copy of CML, and its spatial 
average w{t) = j;'^f^iw{i,t). To relate this problem 
with BKPZ one argues that in the continuous limit and 
close to the synchronized state, the evolution equation 
for w{i,t) is given as a Langevin equation with multi- 
plicative Gaussian noise which then, using the Hopf-Cole 
transformation h = — ln(u;), is transformed into BKPZ. 
Let us notice that CML's are deterministic systems and 
the noise has only effective meaning, mimicking their 
chaotic behaviour. The problem of a relation of such de- 
terministic systems with stochastic counterparts is very 
interesting and recently is drawing some attention [lOj. 
In the above representation the desynchronized phase in 
CML {w{t) > 0) corresponds to the interface pinned rel- 
atively close to the wall {{hi) < oo). In the synchronized 
phase {w{t) 0) the interface depins and drifts away 
{{hi) oo). Perfectly synchronized state {{w{t)) — 0) is 
reached only after infinitely long time. The above analy- 
sis requires the differentiability of the local map and thus 
is not applicable to discontinuous maps. However, nu- 
merical results show that the relation with BKPZ breaks 
down also for continuous but sufficiently steep maps |^. 
From a theoretical point of view it would be desirable to 
understand why such properties of the local map affect 
the nature of the ST and move into the DP universal- 
ity class. Let us notice, that if the coupled CML system 
enters the synchronized state, it will remain in this state 
forever. Such a state can be thus considered as an absorb- 
ing state of the dynamics, although it cannot be reached 
in any finite time. Well-developed techniques are avail- 
able to study phase transitions in models with absorbing 
states (Tlj| . 

It is clear that the problem of synchronization in ex- 
tended systems is related with a number of very inter- 
esting problems in nonequilibrium statistical mechanics 
such as the KPZ model, nonequilibrium wetting f^l or 
directed percolation. Recently, some arguments were 
given that a notoriously difficult particle system, the so- 
called PCPD model, might be also related with these 
models 13|. Deeper understanding of these problems and 
their mutual relations would be certainly desirable. 

In the present paper we introduce an interfacial, dis- 
crete (SOS) model which is inspired by the dynamics of 
CML with discontinuous maps. The synchronization of 
CML is thus related with a depinning of interface from a 
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binding wall. Numerical calculations show that the uni- 
versality class of the depinning transition in our model 
depends on the choice of binding, which enters the dy- 
namics through certain transition rates. For an exponen- 
tially decaying binding, the depinning belongs to the DP 
universality class. In this case the overall behaviour of 
the model is very similar to Sneppen's model of interface 
propagation in a random environment that is driven by a 
certain extremal dynamics 0, . For a power-law de- 
caying binding the depinning transition is characterized 
by a different set of critical exponents. In the case of a 
rapid decay these exponents are very close to those ob- 
tained for the bosonic version of PCPD model In the 
case of a slow decay, in the entire unbounded phase the 
interface remains critical, that in itself is an interesting 
property of a nonequilibrium system. 

In our model discrete site variables hi = 1,2,... are de- 
fined on a one-dimensional lattice of size L {i = 1, . . . , L) 
with periodic boundary conditions — hi). In an 

elementary update we select randomly a site i and change 
the variable hi, and possibly neighbouring ones, accord- 
ing to the following rule: (i) with probability p{hi) one 
sets hi — hi^i = hi^i — 1; (ii) with probability 1 — p(hi) 
the site variable hi increases by unity (hi — > hi + 1). Dur- 
ing a unit of time L elementary updates are performed. 
To complete the definition we have to specify the function 
p{h). To allow a drift toward h = oo the function p{h) 
must decay to zero for h ^ oo. Numerical results that 
we present below are obtained for two cases: p{h) = e~''^ 
(model I) and p{h) = a{h+l)~^ (model II), where 7 > 
and a > are control parameters of the model. 

To make a link with synchronization in CML's the fol- 
lowing remarks are in order. Numerical calculations for 
CML's with discontinuous maps show that the Lyapunov 
exponent that governs the evolution of the synchroniza- 
tion error w{i,t) is negative in the vicinity of the transi- 
tion 0. Approximately, the evolution of w{i,t) is thus 
made of consecutive contractions {w{i,t) — > cw{i,t) and 
c < 1) that, due to discontinuity of the map, are from 
time to time interrupted by discontinuous changes that 
might substantially increase the value of w{i,t). To no- 
tice a link between the dynamics of synchronization er- 
ror in CML and our model we introduce new variables 
Wi = e"''' (let us notice a similarity to the inverse Hopf- 
Cole transformation). Indeed, the increase of hi by unity 
according to the rule (ii) decreases Wi by a factor e that 
corresponds to the contraction of w{i,t). The first rule 
mimics the discontinuous jumps of w{i,t). Since local 
maps in CML's are coupled, a jump at a site i also af- 
fects its neighbours. 

Monte Carlo simulations of our model are similar to 
those of other models with absorbing states j^J- For 
some details related to the fact that the model needs an 
infinite time to reach an absorbing state see e.g., 0. We 
observed that for sufficiently large 7 the interface depins 
from the h = 1 wall and drifts away. For smaller 7 the 
model remains in the active phase with the interface rel- 
atively close to the wall. To examine the nature of the 
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FIG. 1: The steady state activity lo as a function of 7 for 
model I. The inset shows the logarithmic scaling of w at crit- 
icality. The slope of the dotted line corresponds to the DP 
value P = 0.2765. The simulation time was tsim = 10"" plus 
2 ■ 10* discarded for relaxation. 



phase transition we introduced wit) = (x Si^i ^)) 
that in the steady state is denoted as w (in the following 
we refer to this quantity as activity). Of course in the 
absorbing phase w = and in the active phase w > 0. 
Upon approaching the critical point w typically exhibits 
a power-law decay w ^ — jc)~^ : with a character- 
istic exponent P and for the critical point located at 
7 = 7c. Moreover, we studied the time dependence of 
'w{t). One expects that at criticality this quantity has a 
power law decay w{t) ~ , where O is another char- 
acteristic exponent. We also used the so-called spread- 
ing technique First, we set hi = oo for every but 
one site (io) that was set to unity. Then we monitored 
the subsequent evolution of the model (actually, of in- 
terest are only sites with positive w{i,t) i.e., with finite 
hi) measuring the average activity of the system w(t), 
the survival probability P{t), and the averaged spread 
square R^{t) = -^^^ J2i ^(*' 0(* ~ *o)^- One expects that 
at criticality: w{t) - P{t) - f-^ , and R'^(t) - t^, 
that at the same time defines the critical exponents i], S, 
and z. 

First, let us describe the results obtained for model I 
{p{h) = e~'^'^). From the steady state measurements of 
w (Fig. HJ we estimate 7^ = 0.791(1), and P = 0.28(1). 
Such a location of the critical point is confirmed from 
time-dependent simulations (Fig. |3^) and the spreading 
method (Fig. ^p). From these data we estimate 9 = 
0.160(2), r] = 0.317(5), S = 0.16(1), and z = 1.26(1). The 
results for P{t) and R^{t) are not presented. Obtained 
values of critical exponents clearly show that model I 
belongs to the DP universality class for which ^: P — 
0.2765, r] = 0.3137, = 6 = 0.1595 and z = 1.265. 

To have a more complete insight into the critical be- 
haviour of our model we calculate the average interfacial 
width iy(i) = {[Y.f=i{Hht)-{h{i,t)})rV^'^- At critical- 
ity this quantity typically behaves as W{t) ~ where 
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FIG. 2: (a) The time-dependent activity w{t) as a function of 
time t for (from top) 7 = 0.788, 0.790, 0.791 (critical point), 
0.792, and 0.794 (Model I). The results are averaged over 100 
independent runs and simulations were done for L = 5 ■ 10*. 
(b) The time-dependent activity w{t) as a function of time 
t calculated using the spreading method for (from top) 7 = 
0.787, 0.789, 0.791 (critical point), 0.793, and 0.795 (Model 1). 
The results are averaged over 10^ independent runs. The slope 
of the dotted hne corresponds to the DP value i] — 0.3137 
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FIG. 4: The steady state activity 10 as a function of 7 for 
model II for a — 1. The inset shows the logarithmic scaling 
of w at criticahty. The slope of the dotted line corresponds 
to (3 — 0.36. The simulation time was tnim = 10^ plus 2 • 10* 
discarded for relaxation. 
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FIG. 3: The interface profile for model I at criticahty (7 = 
7c = 0.791) after t = 10* MC steps. The inset shows that the 
interface width W{t) grows linearly in time. 



/?' is the growth exponent. For model I simulations show 
almost linear increase of W{t) with time (inset in Fig.|3I) 
and we estimate /?' = 1.0(1). Such a value of the growth 
exponent indicates very strong fluctuations, which is con- 
firmed through a visual inspection of the interface profile 
(Fig.|2Jl. In principle, from direct calculations for CML's 
we can obtain an interface profile as a logarithm of the 
synchronization error. Our preliminary calculations for 
CML model examined by Ahlers and Pikovsky t] show 
that large fluctuations, similar to those in model I are 
clearly seen. However, to quantify these fluctuations and 
calculate e.g., the width W{t), simulations on a longer 
time scale are required. Since for CML dynamics is de- 
fined in terms of w's rather than /I's, simulations for 
longer time, which must probe states of extremely small 
synchronization error, severely suffer from the finite accu- 
racy of numerical computations. On the other hand our 
SOS model easily allows us to examine such a long-time 



regime. 

It is already known that some SOS-like models bel ong 
to the DP universality class. In particular, Alon et al. [l3| 
introduced a model of interface roughening with edge 
evaporation. In this model the growth of the interface 
width is logarithmic in time and is thus much different 
than in our model. As far as the critical behaviour is 
concerned, our model seems to be more closely related 
with the Sneppen's model of interface spreading in envi- 
ronment with quenched disorder |14| . His model is driven 
by extremal dynamics and has the growth exponent /?' 
close to unity. Later it was established that this model 
is actually equivalent to the DP model sitting exactly at 
the critical point [T5|. Important ingredients of Snep- 
pen's model are quenched disorder and extremal dynam- 
ics. The fact that our model, that misses these features, 
exhibits essentially the same behaviour is in our opinion 
quite interesting and worth further examination. 

Using the same procedure we studied the model II 
{P{h) = a{h + ly). First we kept a = 1 fixed and 
varied only the parameter 7. Some of our numerical re- 
sults are shown in Fig. 21-Fig- IHl Using these data we 
estimate: 7^ = 1.74925(5), /3 = 0.36(1), 9 = 0.185(3), 
ri = -0.03(1), 6 = 0.445(5), and z = 1.19(1). All expo- 
nents considerably differ from DP exponents. To support 
these estimations, let us notice that the hyperscaling re- 
lation I = + (5-f'77is satisfied by the above values. 
It is rather surprising for us to observe that our values 
of exponents /3, O and z are in a very good agreement 
with recent estimation for the so-called bosonic version 
of a Pair Contact Process with Diffusion (PCPD) 
It would be interesting to examine whether this is only 
a numerical coincidence or there is a deeper relation be- 
tween these two problems. 

We also studied model II for fixed 7 and varying the 
amplitude a. Monte Carlo simulations show that for 
7 = 2 the critical behaviour of the model seems to be 
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FIG. 5: (a) The time-dependent activity w{t) as a function of 
time t for (from top) 7 = 1.746, 1.748, 1.74925 (critical point), 
1.750, and 1.752 (Model II, a = 1). The results are averaged 
over 100 independent runs and simulations were done for L — 
5 • 10*. (b) The time-dependent activity U'(t) as a function of 
time t calculated using the spreading method for (from top) 
7 = 1.7488, 1.749, 1.7492, 1.74925 (critical point), 1.7493, 
1.7494, 1.7496 and 1.75 (Model II, a = 1). The resuhs are 
averaged over 10^ independent runs. 




the same as for 7 = 7c = 1.74925 and a = I. It suggests 
that for a certain range of 7 the depinning transition be- 
longs to the same universality class. However, a different 
behaviour was observed for 7 = 1. Indeed, in this case 
our simulations suggest (Fig.lSJ that the depinning tran- 
sition is characterized by the exponents /? ~ 3.0(3) and 
Q = 0.62(5). What is also interesting, the entire un- 
bounded phase a < ttc ~ 0.07 seems to be critical with 
a power- law decaying order parameter w ^ , where 
X ~ 0.8(5), and x might slightly vary with a. Such a 
novel behaviour is also worth further studies. 



One of the future problems would be to check whether 
there are some other types of critical behaviour for this 
kind of SOS models. Clearly, an important ingredient 
that determines the critical behaviour is the form of the 
function p{h). It would be interesting to explore some 



other functions (e.g. 



or [ln(ft,)] for a possibly 



new behaviour. Although we were not able to recover 
the BKPZ universality class within our approach, there 
is still a possibility that for a certain choice oip{h) such a 
behaviour might appear. It would be also interesting to 
check whether a new critical behaviour found in model 
II has a counterpart in a synchronization transition in 
CML. 



FIG. 6; The time-dependent activity w{t) as a function of 
time t for (from top) a = 0.11, 0.09, 0.07 (critical point), 0.05, 
and 0.04 (Model II, 7=1). The results are averaged over 200 
independent runs and simulations were done for L = 5 ■ 10*. 
Inset shows the logarithmic scaling of the the steady-state 
activity ui as a function of a — Oc with — 0.07 calculated 
for L = 2 • 10*(+) and i = 5 • 10*(x) (Model II, 7 = 1). 
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